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Abstract 

In this paper measurements are presented of 7l ± , p and p production at mid-rapidity 
(\y\ < 0.5), in Pb-Pb collisions at y/s^m = 2.76 TeV as a function of centrality. The mea- 
surement covers the transverse momentum (pj) range from 100, 200, 300 MeV/c up to 
3, 3, 4.6 GeV/c, for %, K, and p respectively. The measured pj distributions and yields 
are compared to expectations based on hydrodynamic, thermal and recombination models. 



The spectral shapes of central collisions show a stronger radial flow than measured at lower 
energies, which can be described in hydrodynamic models. In peripheral collisions, the pi 
distributions are not well reproduced by hydrodynamic models. Ratios of integrated parti- 
cle yields are found to be nearly independent of centrality. The yield of protons normalized 
^ to pions is a factor ~ 1 .5 lower than the expectation from thermal models. 
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1 Introduction 

The ultimate goal of heavy-ion collisions is the study of the properties of a deconfined and 
chirally-restored state of matter: the Quark-Gluon Plasma (QGP). Indications of its existence 
have been already provided by previous studies at the SPS [1J and at RHIC |2]-[6]|. With the 
advent of the LHC a new energy regime is being studied, aiming at a precise characterization 
of the QGP properties. 

The matter created in heavy-ion collisions exhibits strong collectivity, behaving as a nearly- 
perfect liquid as observed at RHIC B71 0. Its collective properties can be studied through 
transverse momentum (px) distributions of identified particles. A solid understanding of the 
bulk properties of the expanding fireball is necessary for the interpretation of many observables. 
Any signal produced in the QGP phase has to be folded with the space-time evolution of the 
whole system, which has to be taken into account for comparison of theory and data. 

Models based on relativistic hydrodynamics have been very successful in describing observ- 
ables such as the transverse momentum distributions of identified particles, up to a few GeV/c. 
These distributions contain information about the transverse expansion and the temperature at 
the moment when the hadrons decouple from the system [HlQil. It is commonly assumed that 
a significant fraction of the collective flow builds up in the expansion of the fireball in the ini- 
tial partonic phase 10. In this picture, the system would cool down as a consequence of the 
expansion and undergo a phase transition from a partonic to a hadronic phase. The hadrons 
continue to interact, building up additional collective flow and potentially changing the relative 
abundances. The hadronic yields are fixed at the moment when inelastic collisions no longer 
play a role in the system lflTT - fl4l . However, it is usually assumed that the hadronic phase 
does not affect particle abundances. lfT0l[T5l . It was also suggested that the temperature of the 
hadronic ("chemical") freeze-out can be related to the phase transition temperature [fT2l[T6l[T7l . 
Abundances of particles have been fitted very successfully over a wide range of energies (from 
a/snn = 2 GeV to v/Snn = 200 GeV [HUES EES) with thermal (or "statistical hadronization") 
models. From these fits, one can extract the thermal properties of the system at the moment 
when the particle abundances are fixed, the key parameters being the "chemical freeze-out" 
temperature T c ^ and the baryochemical potential \Xb (determined by the net baryon content of 
the system). As will be discussed, the new data presented in this work seem to question part 
of the assumptions in these models, as also reported in IfTBl . The system eventually decouples 
when elastic interactions cease, at the "kinetic freeze-out" temperature T^ n . This temperature, 
together with the expansion velocity at the moment of decoupling can be inferred from the pj 
distributions of identified particles. 

In the intermediate pi region (2 < pi < 8 GeV/c) the baryon-to-meson ratios have been shown 
to reach values > 1 for pi ~ 3 GeV/c, much larger than in pp collisions |fT9ll . It was suggested 
that this could be a consequence of hadronization via recombination of quarks from the plasma 
(in the coalescence models [|20ll2T1l ). 

In this paper, we present the measurement of pi spectra of 7l ± , K , p and p in Pb-Pb collisions 
at v / snn = 2.76 TeV, as a function of centrality and over a wide pj range (from 100, 200, 300 
MeV/c up to 3, 3, 4.6 GeV/c, for n, K, and p respectively) . The ALICE experiment, thanks 
to its unique Particle IDentification (PID) capabilities, is well suited for these measurements. 
Previous results on identified particle production in pp collisions have been reported in ([22|. 
The paper is organized as follows. In Sec. [2] the data sample and the analysis technique are 
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Fig. 1: ITS standalone (left) and global (right) tracking efficiency for pions, kaons and protons (top to 
bottom) in central and peripheral collisions. 



discussed. The systematic uncertainties are presented in Sec. [3} and the results in Sec.|4} These 
are discussed in the context of theoretical models in Sec. [5] Finally, we come to our conclusions 
in Sec. [61 



2 Data Sample and Analysis Method 

The data used for this analysis were collected during the first Pb-Pb run at the Large Hadron 
Collider in the Autumn of 2010. The sample consists of 4 million events, after event selection. 

A detailed description of the ALICE detector can be found in Il23l 1241 . The central tracking 
and PID detectors used in this work cover the pseudorapidity range \t} \ < 0.9 and are, from the 
innermost one outwards, the Inner Tracking System (ITS), the Time Projection Chamber (TPC), 
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Fig. 2: (color online) Performance of the PID detectors: (top) dE/dx distribution measured in the ITS, the 
continuous curves represent the Bethe-Bloch parametrization, the dashed curves the asymmetric bands 
used in the PID procedure; (middle) dE/dx measured in the TPC with global tracks (see text for the 
definition of global tracks), the continuous curves represent the Bethe-Bloch parametrization; (bottom) 
fit of the TOF time distribution with the expected contributions for negative tracks and for the kaon mass 
hypothesis, in the bin 2.4 < pj < 2.5 GeV/c. 

the Transition Radiation Detector (TRD), and the Time-Of-Flight array (TOF). The detector 
features a small material budget (~ 0.1 X/Xq for particles going through the TPC) and a low 
magnetic field, which allow for the reconstruction of low pj particles. The central detectors 
are embedded in a 0.5 T solenoidal field, whose polarity was reversed to allow for systematic 
studies. 

A pair of forward scintillator hodoscopes, the VZERO detectors (2.8 < rj < 5.1 and —3.7 < 
77 < —1.7), measures the arrival time of particles with a resolution of one ns and was used for 
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triggering purposes and centrality determination. The data were collected using a minimum 
bias trigger requiring a combination of hits in the two innermost layers of the ITS (Silicon Pixel 
Detector, SPD, see below) and in the VZERO. The trigger condition used during the data-taking 
has been changed as a function of time, to cope with the increasing luminosity delivered by the 
LHC. This time dependence was eliminated off-line by requiring two hits in the SPD and one 
hit in either of the VZERO detectors. This condition was shown to be fully efficient for the 90% 
most central events G51l . 

The signal in the VZERO was required to lie in a narrow time window of about 30 ns around the 
nominal collision time, in order to reject any contamination from beam-induced backgrounds. 
Only events produced in the vertex fiducial region of |V Z | < 10 cm were considered in the 
analysis (where V z is the vertex position along the beam direction). In the sample used for this 
analysis, a non-negligible fraction of the ions were located outside of their nominal RF bucket in 
the bunch, giving rise to "satellite" collisions. These events are produced well outside the vertex 
fiducial region, but could give rise to "fake" vertices due to the combinatorial algorithm which 
reconstructs the vertices assuming that particles are coming from the area around the nominal 
region. These events were rejected cutting on the correlation of arrival times of beam fragments 
to a pair of "Zero Degree Calorimeters" (ZDCs), placed close to the beam pipe, 114 m away 
from the interaction point on either side of the detector. For details, see Il26ll27ll . 

The VZERO amplitude distribution was also used to determine the centrality of the events. 
In a first step, it was fitted with a Glauber Monte Carlo model to compute the fraction of the 
hadronic cross section corresponding to any given range of the VZERO amplitude. Based on 
these studies, the data were divided in several centrality percentiles, selecting on signal am- 
plitudes measured in the VZERO l|25ll . The results in this paper are reported in ten centrality 
bins, ranging from 0-5% to 80-90%, The centrality intervals and the corresponding charged 
particle multiplicity measured in |tj| < 0.5 (called dN c h/drf in the following) are summarized 
in sec. |4j The dN c h/drj in the centrality bin 80-90%, shown here for the first time, was com- 
puted with the same analysis as described in Il26ll . The contamination from electromagnetic 
processes is negligible down to 80% centrality. In the bin 80-90% an upper limit for this con- 
tribution was estimated as 6%, using the energy distributions of the ZDCs and looking for the 
single (or few) neutron peaks on top of the distribution which would be expected for hadronic 
interactions 11251 1281. 

The production of n ± , K , p and p was measured at mid-rapidity (|y| < 0.5) via three indepen- 
dent analyses, each one focusing on a sub-range of the total pj distribution, with emphasis on 
the individual detectors and specific techniques to optimize the signal extraction. The ranges 
covered by the three analyses are summarized in Table [T] 

The ITS is composed of six layers of silicon detectors using three different technologies. The 
two innermost layers, based on a silicon pixel technology (SPD), are also used in the trigger 
logic, as they provide online the number of pixel chips hit by the produced particles, as men- 
tioned above. The four outer layers, made of drift (SDD) and strip (SSD) detectors, provide 
identification via the specific energy loss. Moreover, using the ITS as a standalone tracker 
enables the reconstruction and identification of particles that do not reach the TPC (at low mo- 
mentum) or cross its dead sectors. This makes the identification of n, K, and p possible down to 
respectively 0.1, 0.2, 0.3 GeV/c in pj. In the first analysis, "ITS standalone" tracks and dE/dx 
were used. At least 4 points are required to form a track, out of which at least 1 must be in the 
SPD and 3 in the drift or strip detectors. With such a small number of tracking points and the 
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high multiplicity in central heavy ion collisions, the probability of having tracks with wrongly 
associated clusters is not negligible. This contribution is strongly suppressed by applying a cut 
on the % 2 P er cluster < 2.5 and not allowing tracks to share clusters. These cuts, however, intro- 
duce a strong centrality dependence of the efficiency, as shown in Fig. [T] (left). The efficiency 
saturates at ~ 0.6 mostly due to the requirement of an SPD point: this detector was only ~ 80% 
operational in 2010 11291 . For each track, the AE /da is calculated using a truncated mean: the 
average of the lowest two points in case four points are measured, or a weighted sum of the low- 
est (weight 1) and the second lowest point (weight 1/2), in case only three points are measured. 
The final AE /da resolution is about 10%. The particle identity was assigned according to the 
distance, measured in units of the resolution, from the expected energy loss curves. While no 
upper limit on this distance was in general used, a 2o lower bound was applied in the case of 
pions, to remove contamination from electrons at low pj. The procedure results in asymmetric 
ranges around the curves for it, K, and p, to reflect the asymmetric nature of the energy loss 
(Fig. [2], top). The range of this analysis is determined at low pj by the ITS standalone tracking 
efficiency, at high pi by the contamination from other particle species: the analysis is stopped 
when the systematic uncertainty coming from this is no longer negligible. 

The other two analyses were based on global tracks, which combine the information from the 
ITS, the TPC and the TRD. This provides good resolution in the transverse distance of closest 
approach to the vertex, DCA X> , and hence good separation of primary and secondary particles. 
The track selection required at least 70 clusters in the TPC and at least 2 points in the ITS, out 
of which at least one must be from the SPD, to improve the DCA Yy resolution. The tracking 
efficiency, shown in Fig. [T] (right), depends only mildly on centrality. It saturates at ~ 70% 
because of the inactive channels in the SPD. Simulations with a fully operational ITS show that 
the intrinsic efficiency of the detector is > 90%. The rise of the efficiency at low pj is due to 
interactions with the detector material, and it is thus sharper for protons. The efficiency reaches 
a maximum when the curvature is big enough so that a track can cross the TPC readout chamber 
boundaries within a relatively small area so that the two track parts can be easily connected (at 
pi ~ 0.6 GeV/c). The straighter tracks at higher pj are affected by the geometrical acceptance. 
This effect is more pronounced for protons than for pions, because the efficiency is folded with 
the decay probability for pions. While protons are stable particles, there is a non-zero decay 
probability for pions which decreases towards higher pj. The shape of the efficiency for kaons, 
on the other hand, is dominated by the decay probability. 

The TPC identifies particles via the specific energy loss in the gas: up to 159 samples can be 
measured. A truncated mean, utilizing only 60% of the available samples, is employed in the 
analysis (Fig. [2j middle). This leads to a Gaussian (and hence symmetric) response function, 
in contrast to the ITS. The resolution is ~ 5% in peripheral and ~ 6.5% in central collisions. 
Further outwards, the TOF measures the time-of-flight of the particles, allowing identification 
at higher pj. The TRD tracking information, if present, is used to constrain the extrapolation 
to the TOF. The total time resolution is about 85 ps and it is determined by three contributions: 
the intrinsic timing resolution of the detector and associated electronics, the tracking and the 
start time. This makes the identification possible out to pi = 3 GeV/c for pions and kaons, and 
4.6 GeV/c for protons. 

In the intermediate pj range, track-by-track identification is possible, based on the combined 
TPC and TOF signals ("TPC/TOF"). It was required that the particles are within 3a of the 
expected values measured in the TPC and/or TOF. The TOF information was used starting at 
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Analysis 




K 


P 


ITS standalone 


0.10-0.60 


0.20-0.50 


0.30-0.60 


TPC/TOF 


0.20-1.20 


0.25-1.20 


0.45-1.80 


TOF fits 


0.50-3.00 


0.45-3.00 


0.50-4.60 



Table 1: pj range (GeV/c) covered by the different analyses. 





-3-2-10 1 2 3 



DCA xy (cm) 

Fig. 3: (color online) DCA V> of protons (top) and antiprotons (bottom) in the /jj-range between 0.6 
GeV/c and 0.65 GeV/c together with the Monte Carlo templates which are fitted to the data (0-5% 
most central collisions). The dashed areas represent the individual templates, the continuous curve the 
combined fit. 
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Pt = 0.65, 0.6, 0.8 GeV/c for n, K, and p respectively, where the 3(7 compatibility cut was 
required for both TPC and TOR The additional requirement of a matching TOF hit reduces the 
overall efficiency shown in Fig. [T] (right) by about 30%, due to the TOF geometrical acceptance 
and to the additional material (see also Sec. [3]). The range of this analysis is determined at low 
Pt by the global tracking efficiency, at high pi by the contamination from other particle species. 
Finally, in the third analysis, a statistical identification based on the TOF signal alone was used 
("TOF fits"). In order to extend the measurement beyond the region of clear separation, the 
difference between the measured time-of-flight and the expected one for the particle species 
under study (normalized to the resolution) was examined. For each pj bin, this distribution was 
fitted with the expected shapes (called "templates" in the following) for n, K, and p, allowing the 
three particles to be distinguished when the separation is as low as ~2c An additional template 
is needed to account for the tracks which are associated to a wrong hit in the TOF (mismatch). 
The templates are built from data, using the measured TOF time response function (described 
by a Gaussian with an exponential tail) and sampling real tracks to get a realistic track-length 
distribution. This fit was repeated for each mass hypothesis, to allow for the calculation of 
the correct rapidity interval. A worst-case example, for 0-5% central collisions, in the bin 
2.4 < pi < 2.5 GeV/c and for the kaon mass hypothesis, can be seen in Fig. [2] (bottom). 

In this paper, results for "primary" particles are presented, defined as prompt particles pro- 
duced in the collision, including decay products, except those from weak decays of strange 
particles. The contamination from secondary particles produced by weak decays or interac- 
tion with the material is mostly relevant for (anti)protons. Since strangeness production is 
typically underestimated in current event generators, and low pi interactions with the material 
are not modeled perfectly in transport codes, the contamination was extracted from data. The 
transverse distance-of-closest-approach (DCAxy) distribution for selected tracks was fitted with 
three distributions ("templates" in the following) corresponding to the expected shapes for pri- 
mary particles, secondaries from material and secondaries from weak decays, as extracted from 
Monte Carlo. This contribution can reach 35% at pi = 300 MeV/c for protons. An example 
of DCA xy fits for the bin 0.6 < pr < 0.65 GeV/c is shown in Fig. [3j the shapes of the three 
contributions are very different. The distribution corresponding to primary particles reflects the 
resolution of the DCA xy . The small non-Gaussian tails seen in the figure are mostly due to large 
angle scattering [|30l . tracks with wrongly assigned ITS clusters or to the combination of tracks 
with hits on the first and/or second ITS layer which have slightly different DCA xy resolution. 
The secondary particles from weak decays show a wider distribution which reflects the large cx 
of weakly decaying particles (of the order of several centimeters). Finally, the secondaries from 
material show very flat tails at high DCAxy. The last contribution is negligible in the case of 
antiprotons. The DCAxy distribution for pions is similar to the one for protons, but with much 
suppressed contribution of secondaries. The distribution for kaons is almost entirely composed 
of primary particles. 

The efficiency correction and the templates used in the secondary correction procedure were 
computed with about 1 million Monte Carlo events, generated using the HIJING OTI event gen- 
erator, tuned to reproduce approximately the dN c h/drf as measured for central collisions [|27l . 
The transport of particles through the detector was simulated using GEANT 3 11321 . The results 
of the three analyses were combined using the (largely independent) systematic uncertainties as 
weights in the overlapping ranges, after checking for their compatibility. 
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3 Systematic Uncertainties 

Several effects related to the transport of particles through the detector and support materials 
contribute to the systematic uncertainties, namely: (i) the amount of secondary protons which 
are produced in the material, (ii) corrections for energy loss through electromagnetic interac- 
tions, (iii) absorption of antimatter due to hadronic interactions. The first effect is taken into 
account in the data-driven DCA xy fit procedure, and its uncertainties are estimated with the 
DCAxy fit variations discussed below. The uncertainties due to the second effect were estimated 
varying the material budget in the simulation by ± 7% and are of the order of 3% and 1%, 
for pions and kaons, 5% and 2% for protons in the two lowest px-bins respectively and then 
quickly decrease further, becoming negligible towards higher momenta. In order to account for 
the last effect, different transport codes (GEANT 3, GEANT 4 [33] and FLUKA fl34l) were 
compared. GEANT 3 (the default transport code used in this work) is known not to reproduce 
the cross sections relevant for the interactions with the material at low pi 051 . The efficiencies 
were scaled with a factor computed with a dedicated FLUKA 0511 simulation. This allows the 
uncertainty for anti-protons to be limited to < 6% despite the large absorption cross-section, 
as there are sufficient existing data on p— nucleus collisions to validate the transport codes. A 
detailed comparison of GEANT 3 with the few existing measurements of hadronic interaction 
cross sections of low momentum kaons and pions O6T4401 reveals differences of about 20-30%. 
After folding with the relevant percentage of particles which are lost due to hadronic interaction 
before entering the TPC or the TOF, the resulting uncertainty is of the order of 2-3% for K 
and % and below 1% for K + . 

Uncertainties in the estimate of contamination from secondary particles can arise from dif- 
ferences in the DCA Yy distributions between data and Monte Carlo. In particular, the DCA Yy 
distribution for secondary particles coming from weak decays is affected by the ct of the decay- 
ing mother, and the actual template used in the fits is a mixture of contributions from different 
particles. The uncertainties due to the secondary subtraction procedure were estimated for all 
analyses by varying the range of the DCAxy fit, by using different track selections (for instance 
using TPC-only tracks), by applying different cuts on the (longitudinal) DCA Z and by varying 
the particle composition of the Monte Carlo templates used in the fit. Overall, the effect is of the 
order of 3% to 1% for pions and 4% to 1% for protons. Such an uncertainty is not relevant for 
kaons, which have a negligible contamination from secondary particles. Moreover the agree- 
ment between the ITS standalone and the TPC analyses provides an additional crosscheck of 
this procedure, because different tracking methods have different sensitivity to contamination 
from secondaries. 

At the lowest pj, the main contribution to the systematic uncertainties comes from the ITS stan- 
dalone tracking efficiency, due to the small number of tracking points and the strong centrality 
dependence. The uncertainty was estimated using the global tracking as a reference. For each 
reconstructed global track, a corresponding ITS standalone track is sought in a narrow window 
defined as At] < 0.03, A0 < 0.03, and Apj < 0.1 x p T . If a track is found, the ITS standalone 
tracking is considered efficient. This pseudo-efficiency is then compared in data and Monte 
Carlo to extract the corresponding systematic uncertainty, estimated to be about 10% for central 
and 3% for peripheral collisions. The efficiency of the ITS standalone reconstruction depends 
strongly on centrality (Fig. [T] left), so that a careful matching between the multiplicity in data 
and Monte Carlo is important. As mentioned, the simulation was adjusted to reproduce the 
d/V ch /dT] measured for central collisions. The occupancy of all ITS layers and the reconstructed 



9 



Centrality Dependence of %, K, p in Pb-Pb at v /Snn = 2.76 TeV 



ALICE Collaboration 



Table 2: Main sources of systematic uncertainty. 



effect 


n ± 




p and p 


pj range (GeV/c) 


0.1 3 


0.2 3 


0.35 


4.5 


correction for 
secondaries 


1.5% 1% 


negl. 


4% 


1% 


material 
budget 


5% negl. 


3% negl. 


3% 


negl. 


haHronip 
interactions 


2% 1% 


4% 1% 


6% 
4% 


l%(p) 
negl. (p) 


pj range (GeV/c) 


0.1 0.5 


0.2 0.5 


0.35 


0.65 


ITS tracking 

efficiency (Central 0-5%) 


10% 


10% 




10% 


ITS tracking 

efficiency (Peripheral 80-90%) 


4% 


4% 




4% 


ITS PID 


2% 


4% 




4.5% 


pi range (GeV/c) 


U.J U.UJ 


U.J U.O 


0.5 


0.8 


global tracking 
efficiency 


4% 


4% 




4% 


TPC PID (Central, 0-5%) 


3% 


5% 




1.5% 


TPC PID (Peripheral 80-90%) 


1.5% 


3.5% 




1% 


pj range (GeV/c) 


0.5 3 


0.5 3 


0.5 


4.5 


TOF matching 
efficiency 


3% 


6% 




3% 


TOF PID 


2% 7% 


3% 15% 


5% 


25% 



track multiplicity were compared between data and Monte Carlo. Residual differences, due to 
the different dependence of multiplicity on centrality in data and Monte Carlo, contribute to the 
systematic uncertainty on the tracking mentioned above and are of the order 2% . The sensi- 
tivity to the variation of the track selection cuts (DCAxy, % 2 , number of clusters) was found to 
be of the order of 7% for all particles. Moreover, the Lorentz force causes shifts of the cluster 
position in the ITS, pushing the charge in opposite directions when switching the polarity of the 
magnetic field of the experiment (ExB effect). This effect is not fully reproduced in the Monte 
Carlo simulation and was estimated analyzing data samples collected with different magnetic 
field polarities. This uncertainty is only relevant at the lowest pj and is 3% for pions, and 1% 
for kaons and protons. 

The uncertainties related to the ITS PID method were estimated by using different techniques 
for the identification. A 3o cut on the maximum difference in energy loss of the particles was 
applied (instead of using all particles as in the default strategy). Alternatively, an unfolding 
method was used, in which for each pj bin, the dE/dx distribution is fitted with Monte Carlo 
templates for each one of the species, similar to what is done in the TOF fits. The effect is found 
to be 2%, 4%, 7% for n, K, and p respectively. Close to the high pj boundaries of the analyses 
(Table [T]), the possibility of misidentifying other species is not completely negligible. This is 
corrected for with Monte Carlo, and the corresponding uncertainty, due to residual differences 
in particle ratios and dE/dx between data and Monte Carlo is estimated to be 2% for kaons and 
1% for protons. 

The uncertainty in the global tracking efficiency was estimated comparing the track matching 
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Fig. 4: (color online) Transverse momentum (pj) distribution of n, K, and p as a function of centrality, 
for positive (circles) and negative (squares) hadrons. Each panel shows central to peripheral data; spectra 
scaled by factors 2" (peripheral data not scaled). Dashed curves: blast-wave fits to individual particles; 
dotted curves: combined blast-wave fits (see text for details). Statistical (error bars) and systematic 
(boxes) uncertainties plotted. An additional normalization uncertainty (Table 4.2 1 has to be added in 
quadrature. 
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efficiency from TPC to ITS and from ITS to TPC in data and Monte Carlo. The effect is found 
to be ~4%. The uncertainties related to the track selection were investigated by a variation of 
the track cuts. They are estimated (for most central collisions) to be 3% for pions and kaons and 
4% for protons. This uncertainty decreases slightly for peripheral collisions. The uncertainties 
related to the TPC/TOF PID procedure were estimated varying the PID cut between 2 and 5 
o. As pions are the most abundant particles, their corresponding systematic uncertainty does 
not exceed 3% even in regions where the mean pion energy loss crosses the kaon and proton 
curves. On the other hand, it can reach up to 4-5% for kaons. The corresponding systematics of 
the better separated protons are below 1.5%. As the dE/dx resolution becomes slightly worse 
with increasing multiplicity, the systematic uncertainty shows a similar slight increase. The 
values quoted here represent upper bounds for all centralities as they correspond to the 0-5% 
most central collisions. 

The uncertainties on the response functions of the PID detectors were found to have a negligible 
effect: they were carefully tuned using high statistics data and this analysis is limited to regions 
were the separation of different particles is >2a. 

As mentioned above, the tracks reaching the TOF detector have to cross a substantial amount of 
additional material budget (about 23% X/Xq), mostly from the Transition Radiation Detector 
(TRD) [1241 . The systematic uncertainties on the TOF matching were estimated comparing the 
matching efficiency evaluated in Monte Carlo and from data using samples of cleanly identified 
particles in TPC. Good agreement is observed in case of pions and kaons, with deviations at the 
level of at most 3% and 6%, respectively, over the full pi range. 

In the case of protons and antiprotons good agreement is also observed, for pj > 1 GeV/c, with 
maximal deviations of 7.5%. Since the TRD was not fully installed in 2010, the analysis was 
repeated for regions in azimuth with and without installed TRD modules, allowing to cross- 
check the uncertainty on the material. The effect is found to be 3%, 6%, 3% for n, K, and 
p respectively. The uncertainties related to the identification in the "TOF fits" analysis were 
estimated varying the parameters of the expected sources in the fit by ± 10%. This is one of the 
main sources of uncertainty and increases with pj, as the separation between different species 
decreases. It is < 5% at p 7 ~ 0.5 GeV/c and < 7%, 15%, 25% at high p T , for %, K, and p 
respectively. 

The main sources of systematic uncertainties are summarized in Table [2] 

The measured spectra are further affected by an additional normalization uncertainty, coming 
from the centrality estimation. The centrality percentiles are determined with sharp cuts on the 
VZERO amplitude distribution, which are affected by a 1% uncertainty [|25ll . This translates 
into a normalization uncertainty on the spectra. The total normalization uncertainty is about 
+ 12% -8.5% for peripheral events (also including a 6% contribution due to contamination from 



electromagnetic processes) and negligible for central events, see Table 4.2 



4 Results 

4.1 Transverse Momentum Distributions 



The combined spectra in the centrality bins reported in Table 4.2 are shown in Fig. |4} The 



distributions of positive and negative particles are compatible within uncertainties at all pj, as 
expected at LHC energies. This is clearly seen in the ratio of negative to positive pj spectra 
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Fig. 5: (color online) Ratios of negative to positive particles as a function of pj. 



13 



Centrality Dependence of 7t, K, p in Pb-Pb at v /Snn = 2.76 TeV 



ALICE Collaboration 



o 

> 

01 

CD 



I— 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




7f + 71 

0-5% 
123 20-30% 
W\ 70-80% 




0.5 1 1.5 2 2.5 3 

p r (GeV/c) 



> 

CD 

CD 



> 

CD 

CD 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 


3 
2.5 

2 
1.5 

1 

0.5 




0.5 



1.5 



K + + K" 

20-30% 

70-80% 



2.5 



P T (GeV/c) 



P + P 

OH 0-5% 
HI 20-30% 

Em 70-80% 



P T (GeV/c) 



Fig. 6: (color online) Local slopes of the pj distributions of %, K, and p (summed charges) as a function 
of pj and centrality. See text for details. 

shown in Fig. [5] For this reason, we focus in the following mostly on results for combined 
charges. 

The change of shapes with centrality is apparent in Fig.|4j the spectra get harder with increasing 
centrality. The spectra for all particle species have an almost exponential shape at high px m 
central collisions: contrary to what is observed in pp collisions ll22ll there is no obvious onset 
of a high-px power- law tail. In more peripheral collisions, the onset of the power-law is visible 
at high pj. This is quantified in Fig. [6] which shows the local inverse slope 7] oc of the spectra, 
as a function of pj, computed with a fit using five bins in the vicinity of each pj bin with the 
function 



1 dN 
pjdpr 



-Pt/T\ 



(1) 
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This approach guarantees a numerically stable extraction of the local inverse slope parameter 
7] oc , but leads to a significant correlation of the depicted uncertainties for neighboring points. 
Residual correlations in the bin-by-bin systematics artificially increase the uncertainties on 7] oc 
shown in Fig. [6] As seen in the figure, the inverse slopes of K and p become larger with 
decreasing pj (the spectra are flatter at low pj), and this change is more pronounced for central 
collisions. The proton spectrum at low pi is nearly flat, making 7i oc unconstrained (for a flat 
distribution, 7] oc — > °°). Above a certain pj, roughly 1 GeV/c for K and 2 GeV/c for p, the slopes 
do not change with pj for central and semi-central collisions, consistent with an exponential 
shape. Protons and kaons converge to similar values of 7i oc ~ 0.45 GeV/c. For peripheral 
collisions, a modest increase of 7i oc is seen at the highest pj, suggesting the onset of a power- 
law behavior. The 7] oc trend is different for pions. The % spectra are not purely exponential. 
At high pj the power-law rise is more suppressed in central collisions compared to peripheral 
ones. At low pj T\ oc increases with pj, opposite to the trend observed for protons and kaons. 
This steepening of the pion spectra is a general feature of heavy-ion collisions, and is due to the 
large contribution of resonance decays to the pion spectrum, as already noted in [4TJ . Above 
pi ~ 1 GeV/c, the rate of increase of 7i oc is slower than at lower pj, and is less pronounced for 
central collisions. 

4.2 px-Integrated Yields and Mean Transverse Momentum 

In order to extrapolate to zero pjfor the extraction of px-integrated yields and {pr), the spectra 
were fitted individually with a blast- wave function [J9j : 

1 - ^rdrmx/o ( ^^-) K x ( m ? C ™ h P ) , (2) 



Pi dp T Jo \ Tki n J V T t 

where the velocity profile p is described by 

p=tanh- 1 p\ = tanh- 1 ^Q"/3,j . (3) 

Here, mj = ^Jpj + m 2 is the transverse mass, Iq and ^fi the modified Bessel functions, r is the 

radial distance in the transverse plane, R is the radius of the fireball, p\ is the transverse expan- 
sion velocity and j8 s is the transverse expansion velocity at the surface. From these equations 
one can also derive the average transverse expansion velocity (j&r) . The free parameters in the 
fit are the freeze-out temperature T^ n , the average transverse velocity (p\) and the exponent 
of the velocity profile n. This function describes very well all particle species over the whole 
measured pi range (as individual fits). It should be noted, however, that from fits to a single 
particle species no physics meaning can be attached to those parameters. A combined fit to 
different particle species can provide insight on the freeze-out parameters, and this is discussed 
in detail in the next section. The fraction of extrapolated yield is small: about 7%, 6%, 4% for 
K, K, p, respectively. The systematic uncertainties due to the extrapolation amounts to 2.5%, 
3%, 3% for the yields (independent of centrality) and to 2%, 2%, 3% (peripheral) or 2%, 3.5%, 
3.5% (central) for (pr). This was estimated using different fit functions [42J (Boltzmann, mx 
exponential, pj exponential, Tsallis-Levy, Fermi-Dirac, Bose-Einstein), restricting the fit range 
to low pj for those functions not giving a satisfactory description of the spectra over the full 
range. The (px) is computed extrapolating with the blast wave function to 100 GeV/c (infinity, 
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effectively). The difference between the (pj) computed with and without the extrapolation at 
high pi is < 1%, ~ 1.5%, < 1% for n, K, and p respectively. The extracted particle yields and 
(pt) as a function of centrality are summarized in Table 4.2 and Table |4j 



Figure [7] shows the mean transverse momentum (pj) as a function of dN c h/dr}, compared to 
previous results atRH!CQ[|42l|43l|45ll. The (p T ) increases with centrality and is higher than the 
lower energy results for comparable charged particle densities. The (px) at RHIC was found to 
be compatible with a scaling as a function of dA^h/dr] for different energies [42J. This scaling 
is clearly excluded at the LHC. 

The ratios of negative to positive particle yields (Fig. [8]) are compatible with unity for all cen- 
tralities. The p/p ratio, in particular, confirms the expectation of a vanishing baryon transport 
to mid-rapidity at the LHC, in contrast to to the RHIC energy regime, where the p/p ratio was 
found to be about 0.8 at ^/snn = 200 GeV [|42l . The effect of the different antibaryon to baryon 
asymmetry between the two energies is almost absent in the sum of the positive and negative 
charges. Therefore, the ratios K/n = (K + + K - )/ (n + + n~) and p/n = (p + p)/ (n + + n~) 
are compared to RHIC 11421 l44l 1451 in Fig. [9] as a function of charged particle multiplicity. The 
K/ % ratio hints at a small increase with centrality following the trend from lower energy data. 
The p/n ratio suggests a small decrease with centrality and is slightly lower than the RHIC 
measurements. 



5 Discussion 

5.1 Transverse Momentum Distributions and Hydrodynamics 

The change with centrality of the spectral shapes shown in Fig. [4] and Fig. [6] can be interpreted 
in terms of hydrodynamics. A flattening of the spectra, more pronounced at low pi and for 
heavier particles, is expected in the hydrodynamical models as a consequence of the blue-shift 
induced by the collective expansion. The low-px change of the local slope shown in Fig. [6} 
more pronounced for the proton spectra, thus suggests a progressively stronger radial flow with 
increasing centrality. The fact that the inverse slope converges to the same value for p and K at 
high pj is also a generic feature of the blast wave parameterization. 

In order to quantify the freeze-out parameters at ^/snn = 2.76 TeV, a combined fit of the spectra 
with the blast- wave function Eq. ([2]) was performed, in the ranges 0.5-1 GeV/c, 0.2-1.5 GeV/c, 
0.3-3 GeV/c for n, K, p, respectively. The pions at low pj are known to have a large contribu- 
tion from resonance decays, while at high pj a hard contribution (not expected to be described 
by the blast-wave) may set in. Therefore, the values of the parameters extracted from the fit, 
and especially T^ n , are sensitive to the fit range used for the pions. Forcing all species to de- 
couple with the same parameters also makes the interpretation of the results arguable: different 

'The RHIC data are plotted as a function of cW c h/dT] using the measured dN c ii/drj of each individual ex- 
periment [42-44]. In the case of PHENIX, the dN c h/drj are published in 5% percentiles, while the spectra are 
published mostly in 10% percentiles: whenever needed, the value of the dN c h/dr] used in the figure is a linear 
interpolation of the 5% percentiles; the dN c h/dri measurement, moreover, is only available up to 70% centrality 
in PHENIX. Since there is some disagreement in the dN c h/dT] measurements from different RHIC experiments 
for peripheral event, we decided not to plot in Fig. [JJthe PHENIX results below 70% centrality. The discrepancy 
is also visible in Fig. [15] where the difference in normalization between STAR and PHENIX is apparent. The 
STAR (anti)proton measurement is inclusive of products from weak decays of strange particles, and therefore not 
included in all comparisons shown in this section. 



16 



Centrality Dependence of 7t, K, p in Pb-Pb at v /Snn = 2.76 TeV 



ALICE Collaboration 




10 
dA/ ct /dJ7 



> 

CD 

(5 




10 
6NJ6t! 



> 

CD 

(5 



2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 



10 



ALICE, Pb-Pb, v^jm = 2.76 TeV 
STAR, Au-Au, \js NN = 200 GeV 
PHENIX, Au-Au, \!sZ. = 200 GeV 



Full: positive, Open: negative 



P.P 



10 £ 



10 



Fig. 7: (color online) Mean transverse momentum {pj) as a function of d/V c h/dr/ compared to RHIC 
results at ^/snn = 200 GeV B2l |45l . Negative charge results displaced for better readability. Boxes: 
systematic uncertainties. 
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Fig. 8: (color online) Ratios of negative to positive particles as a function of (W ch /dr/, compared to 
previous results at -vAnn = 200 GeV. 
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Fig. 9: (color online) K/n = (K + + Kr))/(7r + + 7z;~) andp/TT 
of dNch/drj, compared to previous results at ^/snn = 200 GeV. 
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Fig. 10: (color online) Left: Results of blast- wave fits, compared to similar fits at RHIC energies [5]; The 
uncertainty contours include the effect of the bin-by-bin systematic uncertainties, the dashed error bars 
represents the full systematic uncertainty (see text for details), the STAR contours include only statistical 
uncertainties. Right: Comparison of fit results for different fit ranges; the error bars include only the 
effect of the bin-by-bin systematics (see text for details). 
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Fig. 11: (color online) Blast-wave parameters (jSp) and Ta n as a function of cWch/dr/, compared to 
previous results at ^Ann = 200 GeV [5 ] (full systematic uncertainties for both experiments). 

particles can in principle decouple at a different time, and hence with a different (j8t) and T^ n , 
from the hadronic medium, due to their different hadronic cross section. These fits by no means 
replace a full hydrodynamical calculation: their usefulness lies in the ability to compare with a 
few simple parameters the measurements at different ^snn- As will be discussed, the param- 
eters extracted from such a combined fit depend on the range used for the different particles. 
Our standard fit ranges were therefore chosen to be similar to the ones used by the STAR col- 
laboration at -y/sNN = 200 GeV/c at the low pi end. The high pj boundaries were extended to 
higher pj as compared to STAR, since at the LHC it is expected that the shapes are dominated 
by collective effects out to higher transverse momenta. The results of the fit are summarized 
in Tableland shown in Fig. 10 (left) and in Fig. 11 The 1-sigma uncertainty ellipses shown 
in the figure reflect the bin-to-bin systematic uncertainties. The uncertainties shown as dashed 
bars and reported in Table [5] also include systematic uncertainties related to the stability of the 
fit: the effect of the variation of the lower fit bound for pions (to test the effect of resonance 
feed-down), the sensitivity to different particle species (i.e. excluding pions or kaons or protons) 
and to fits to the individual analyses. The value of (j&r) extracted from the fit increases with 
centrality, while decreases, similar to what was observed at lower energies (Figs. [10] and [TT 



Tab. [5]). This was interpreted as a possible indication of a more rapid expansion with increasing 
centrality [5J. In peripheral collisions this is consistent with the expectation of a shorter lived 
fireball with stronger radial gradients IfTOl . 

The value of the n parameter, Eq. ([3]), is about 0.7 in central collisions and it increases towards 
peripheral collisions. The large values in peripheral collisions are likely due to the spectrum not 
being thermal over the full range: the n parameter increases to reproduce the power law tail. 

In order to further test the stability of the fit, it was repeated in the ranges 0.7-1.3 GeV/c, 0.5- 
1.5 GeV/c, 1-3 GeV/c ("high p T ") and 0.5-0.8 GeV/c, 0.2-1 GeV/c, 0.3-1.5 GeV/c ("low p r ") 



for 7t, K, p respectively. The effect of the fit ranges is demonstrated in Fig. 10 (right). As can 



be seen, while the value of (jSp) is relatively stable, especially for the most central bins, the 
value of 7ki n is strongly affected by the fit range, with differences of order 15 MeV for the most 
central events. For most peripheral events, (j3t) shows some instability. However, it should be 
noticed that this parameter is mostly fixed by the low-px protons and its uncertainty increases 
significantly when the fit range for the protons is reduced. 

The combined fits in the default range are shown in Fig.[4j compared to fits to individual particle 
spectra. The dotted curves represent the combined blast-wave fits, while the dashed curves are 
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Fig. 12: (color online) Ratio of the measured spectra to the combined blast-wave fit as a function of pj. 
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Fig. 14: (color online) Spectra of particles for summed charges in the centrality bin 20-30% compared 
to hydrodynamical models and results from RHIC at i/5nn = 200 GeV B2l l45ll . See text for details. 
Systematic uncertainties plotted (boxes); statistical uncertainties are smaller than the symbol size. 
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Fig. 15: (color online) Spectra of particles for summed charges in the centrality bin 70-80% compared 
to hydrodynamical models and results from RHIC at i/5nn = 200 GeV ll42l l45l . See text for details. 
Systematic uncertainties plotted (boxes); statistical uncertainties are smaller than the symbol size. 
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the individual fits to each species. The individual fits reproduce the spectra over the full pj 
range in all centrality bins, thus allowing for a reliable estimate of the (pj) and integrated 
yields, as mentioned above. The extrapolation at low and high pj of the combined fits gets 
progressively worse for decreasing centrality. The ratios of the spectra to the combined fits 



are shown in Fig. 12 If the behavior of the spectra is purely hydrodynamical over the full pj 
range considered, one would expect that the fireball blast-wave parameters determined by a fit 
in a limited pj range are able to predict the full shape. This is what is observed in the most 
central bin for protons and kaons. The same is not true for the more peripheral bins, and the pj 
threshold at which the function deviates from the data decreases with centrality, indicating the 
limit of applicability of the hydrodynamical picture (as also discussed below). 

The spectra for summed charges in central (0-5%), semi-central (20-30%) and peripheral (70- 
80%) collisions are compared to hydrodynamical models and previous results in Au-Au colli- 



sions at -y/sNN = 200 GeV in Figs. 13 14 and 15 



A dramatic change in spectral shapes from RHIC to LHC energies is observed, with the protons 
in particular showing a flatter distribution. A comparison between the two energies based on 
the values of (j8t) and T^ n from the combined blast- wave fits Q is shown in Fig. 10 an dJTTj 
For central collisions, about 10% stronger radial flow than at RHIC is observed at the LHCq 



The models shown in Figs. 13 14, and 15 give for central collisions a fair description of the 
data. In the region p T < 3 GeV/c (Krakow @fD), p T < 1.5 GeV/c (HKM02) and p T < 3 GeV/c 
(EPOS ll48l . with the exception of protons which are underestimated by about 30% at low pj), 
the models describe the experimental data within ~20%, supporting a hydrodynamic interpreta- 
tion of the pi spectra in central collisions at the LHC. VISH2+ 1 [|49ll is a viscous hydrodynamic 
model that reproduces fairly well the pion and kaon distributions up to pj ~ 2 GeV/c, but it 
misses the protons, both in shape and absolute abundance in all centrality bins. In this version 
of the model the yields are thermal, with a chemical temperature = 165 MeV, extrapolated 
from lower energies. The difference between VISH2+1 and the data are possibly due to the 
lack of an explicit description of the hadronic phase in the model. This idea is supported by the 
comparison with HKM [|47l|50]|. HKM is a model similar to VISH2+1, in which after the hy- 
drodynamic phase particles are injected into a hadronic cascade model (UrQMD), which further 
transports them until final decoupling. The hadronic phase builds up additional radial flow and 
affects particle ratios due to the hadronic interactions. As can be seen, this model yields a better 
description of the data. The protons at low pj, and hence their total number, are rather well 
reproduced, even if the slope is significantly smaller than in the data. Antibaryon-baryon anni- 
hilation is an important ingredient for the description of particle yields in this model [|47l l50l . 
The Krakow OTl [52l model, on the other hand, uses an ansatz to describe deviation from equi- 
librium due to viscous corrections at freeze-out, which seems successful in reproducing the 
data. A general feature of these models is that, going to more peripheral events, the theoretical 



curves deviate from the data at high pj (Figs. 14 and 15). This is similar to what is observed 



in the comparison to the blast- wave fits, and shows the limits of the hydrodynamical models. 
As speculated in [|46ll . this could indicate the onset of a non-thermal (hard) component, which 
in more peripheral collisions is not dominated by the flow-boosted thermal component. This 

2 The full systematic uncertainties on the parameters, quoted by both STAR and ALICE, include a number of 
checks on the stability of the fit which have a similar effect at different energies. When comparing results of 
different experiments, only the statistical and systematic uncertainties of the pj distributions should be considered, 
but not those related to the stability of the fit. 
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Centrality 


(fr) 


r kin (GeV/c) 


n 


X /ndof 


0-5% 


0.651 ±0.004 ±0.020 


0.095 ±0.004 ±0.010 


0.712±0.019±0.086 


0.15 


5-10% 


0.646 ±0.004 ±0.023 


0.097 ±0.003 ±0.011 


0.723±0.019±0.116 


0.20 


10-20% 


0.639 ±0.004 ±0.022 


0.099 ±0.004 ±0.011 


0.738±0.020±0.118 


0.19 


20-30% 


0.625 ±0.004 ±0.025 


0.101 ±0.004±0.012 


0.779 ±0.022 ±0.133 


0.22 


30^10% 


0.604 ±0.005 ±0.022 


0.106 ±0.004 ±0.012 


0.841 ±0.025 ±0.168 


0.22 


40-50% 


0.574 ±0.005 ±0.016 


0.112±0.004±0.013 


0.944 ±0.029 ±0.142 


0.22 


50-60% 


0.535 ±0.007 ±0.018 


0.1 18 ±0.004 ±0.014 


1.099 ±0.038 ±0.187 


0.28 


60-70% 


0.489 ±0.008 ±0.024 


0.129 ±0.005 ±0.017 


1.292 ±0.052 ±0.194 


0.36 


70-80% 


0.438 ±0.011 ±0.039 


0.1 39 ±0.005 ±0.027 


1.578 ±0.081 ±0.205 


0.40 


80-90% 


0.357±0.016±0.084 


0.151 ±0.006 ±0.044 


2.262±0.191 ±0.498 


0.52 



Table 5: Results of the combined blast-wave fits, in the default ranges 0.5-1 GeV/c, 0.2-1.5 GeV/c, 0.3- 
3 GeV/c for %, K, p respectively. The first uncertainty in the table includes the effect of the bin-by-bin 
systematic uncertainties, the second is the full systematic uncertainty. See text for details. 



picture is further substantiated by the change in the local slopes as seen in Fig. [6| 

The EPOS (2.17v3) model Il48l aims at describing all pj domains with the same dynamical 
picture. In this model, the initial hard scattering creates "flux tubes" which either escape the 
medium and hadronize as jets, or contribute to the bulk matter, described in terms of hydrody- 
namics. After hadronization, particles are transported with a hadronic cascade model (UrQMD). 
EPOS shows a good agreement with the data for central and semi-central collisions. A calcu- 
lation done with the same model, but disabling the late hadronic phase, yields a significantly 
worse description ll48l . indicating the important role of the late hadronic interactions in this 
model. An EPOS calculation for peripheral collisions was not available at the time of writing, 
but it will be important to see how well the peripheral data can be described in this model, since 
it should include all relevant physics processes. Several other models implementing similar 
ideas (hydrodynamics model coupled to a hadronic cascade code, possibly with a description of 
fluctuations in the initial condition) are available in the literature Il53ll54ll but not discussed in 
this paper. The simultaneous description of additional variables, such as the v n azimuthal flow 
coefficients within the same model, will help in differentiating different hydrodynamical model 
scenarios. 



Figure 16 shows the p/n = (p + p)/(7t + + n~) and K/n = (K + +Kr)/(7T+ + n~) ratios as a 
function of pj. Both ratios are seen to increase as a function of centrality at intermediate pj 
with a corresponding depletion at low pj (the pi integrated ratios show little dependence on 
centrality, Fig. [9|. The p/it ratio, in particular, shows a more pronounced increase, reaching 
a value of about 0.9 at pj = 3 GeV/c. This is reminiscent of the increase in the baryon-to- 
meson ratio observed at RHIC in the intermediate pi region |fT9l 15511 . which is suggestive of 
the recombination picture discussed in Sec. [T] It should be noted, however, that a rise of the 
ratio with pi is an intrinsic feature of hydrodynamical models, where it is just due to the mass 
ordering induced by radial flow (heavier particles are pushed to higher pi by the collective 



motion). In Fig. 16 a prediction from two of the hydrodynamical models discussed above 11471 
1521 and a prediction from a recombination model 11201 1561 are shown. As can be seen, the ratio 
for central events is reasonably reproduced by the Krakow hydrodynamical model, while HKM 
only reproduces the data up to pj ~ 1 .5 GeV/c and the recombination prediction is higher than 
the data and predicts a flatter trend in the range 2-3 GeV/c. 

This measurement is currently being extended to higher pj by ALICE using the HMPID detec- 
tor (High Momentum Particle Identification [15711 . a ring imaging Cerenkov), complemented by 
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Fig. 16: (color online) Ratios K/k = (K + + Kr))/(7T + + 7r~) andp/7r = (p+p)/(7T + + 7r~) as a function 
of pj for different centrality bins compared to ratios from the Krakow and HKM hydrodynamic model. 
The p/n is also compared to a recombination model. 
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a statistical identification in the relativistic rise region of the TPC. A complementary study of 
the A/K^ ratio will also provide a good pj coverage and additional constraints. 

5.2 px-Integrated Yields 

The integrated ratios can be interpreted in terms of the thermal models. Figure [9] depicts the 
expectations from these models, which, based on lower energy data, used the values T ~ 
160 - 170 MeV and ji B ~ 1 MeV at the LHC |58l|59]|. Tne K / n ratio is consistent with these 
expectations, while the p/n ratio is found to be lower by a factor of about 1.5. As discussed in 
lTT8ll . this finding was one of the surprises of the first Pb-Pb run at the LHC and still needs to be 
understood. 

Some indication of a similar disagreement between data and the thermal model is also seen 
in the RHIC data, with the proton measurements being 10-20% lower than the thermal model 
predictions |l5l[TT1|. This discrepancy was not considered to be significant, due to experimen- 
tal uncertainties in the subtraction of secondary particles, differences between thermal model 
implementations and model uncertainties Il60"ll . 

5.3 Total proton spectrum 

In order to compare directly with lower energy results, we performed a measurement of protons 
including feed-down from weak decays ("total proton spectrum"), based on tracks reconstructed 
using only TPC information, with no requirements on the ITS. These tracks have a similar effi- 
ciency for both primary and secondary particles, with the difference in the particle composition 
between data and Monte Carlo being of minor importance. The reconstructed sample thus al- 
ready includes most of the secondary particles. The efficiency correction for the total proton 
measurement is about 25% at 450 MeV/c and 5% at 1 GeV/c. In this analysis, the secondaries 
from detector material were subtracted using a DCA xy fit procedure, similar to the one dis- 
cussed in Sec. [2] The total proton spectrum for central collisions (0-5%) is compared in Fig. 17 



with the primary proton spectrum and with the previous STAR measurements of total protons at 
-y/SNN = 200 GeV. The comparison of the spectra confirms the change in shape and yield from 
RHIC to LHC, already discussed for primary protons. The total spectrum was fitted with a blast- 



wave function (also shown in Fig. 17 ) to compute the extrapolation for the extraction of the inte- 
grated yield. Due to the limited p^-coverage of the TPC-only analysis, the fraction of the extrap- 
olated yield amounts to about 25% resulting in a larger extrapolation uncertainty (about 12%) as 
compared to the primary spectrum. We found dN^p/dy = 120 ± 20 (syst + stat) for the com- 
bined spectra of protons and antiprotons, resulting in the ratio p (total)/7T = 0.082 ±0.010. This 
is about 15% lower than the ratio of 0.095 ± 0.01 1 02 measured by STAR at = 200 GeV 
(consistent with what was previously noted for primary protons), but the difference is not sig- 
nificant within uncertainties. 

Subtraction of the total and primary proton yields dN^p/dy = 53 ± 19 (syst + stat) for sec- 
ondary protons. This number, normalized to the primary pion yield, can be compared to the 
thermal model prediction [fT2| with T c h = 164 MeV. There are five main contributions to the 
secondary protons, summarized in Table [6] Rescaling the ratios shown in Table [6]by the mea- 
sured pion yield, we get dVp^ model /dj = 62, in good agreement with the measurement. This 
result suggests that the disagreement between the data and the thermal model is most prominent 
for primary protons, while strange baryons contributing to the production of secondary protons 
are likely better described. If the total protons measured in the data and in the thermal model 
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Particle Decay Channel Branching Ratio p (secondary)/^ 

A p7T 63.9 % 4.42 • 10" 2 

E+ p7T° 51.6% 1.27 • 

S- Atz 99.89 % 5.49 -l(r 3 

E° An 99.53 % 5.58 1(T 3 

£2 AK 67.8 % 9.83 • 10" 4 



Table 6: Contribution to secondary protons estimated in the thermal model 11121 with T c h =164 MeV. 
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Fig. 17: (color online) Transverse momentum pj distribution of total protons based on tracks recon- 
structed using only TPC information (full circles) compared to primary protons (open circles) and cor- 
responding total spectra measured by the STAR collaboration in Au-Au collisions at ^Ann = 200 GeV/c 
(full stars). 0-5% most central events. Box: systematic uncertainty; Statistical uncertainties smaller than 
the symbols; curve: individual blast wave fit. 




are compared, the disagreement gets partially suppressed, because the number of secondaries 
is well reproduced in the thermal model, and the fraction of secondary protons is rather large 
(~ 50%). 

A possible explanation for the difference between the p/n ratio and the predictions from the 
thermal models is antibaryon-baryon annihilation in the hadronic phase P7l I6T1 l62l . The 
p/n ratio depicted in Fig. [9] suggests a decreasing trend with centrality, consistent with the 
antibaryon-baryon annhilation hypothesis: the effect is expected to be less important for the 
more dilute system created in peripheral collisions. It should be noted that all the available 
models incorporating a hadronic phase use the UrQMD If63l [641 hadronic cascade model, and 
the effects induced by the hadronic phase are model dependent. While this microscopic model 
includes annihilation processes, it does not implement the reverse process like nit — > pp. The 
effect of the reverse reactions was investigated in a recent calculation [|65l . It was found to be 
non-negligible, while the net suppression of baryons is still very significant. 

The origin of the low proton yield with respect to the thermal model expectations is not yet 
established and alternative explanations exist in the literature. Implementations of the thermal 
model incorporating non-equilibrium effects predict a reduction of the p/n ratio, although the 
details depend on the exact value of model parameters and actual predictions underpredict the 
measured ratio Il66ll . An alternative explanation may come from the existence of flavor and 
mass dependent pre-hadronic bound states in the QGP phase, as suggested by recent lattice 
QCD calculation and QCD-inspired models [[671 1681 . Additional constraints to discriminate 
between these different scenarios will be provided by a thermal analysis, including in particular 
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strange and multistrange baryons. 

6 Conclusions 

In this paper we presented a comprehensive measurement of %, K, and p production in Pb-Pb 
collisions at = 2.76 TeV at the LHC. Antiparticle over particle ratios are compatible with 
unity at all pj and at all centralities, as expected for LHC. A clear evolution of all spectra with 
centrality is seen, with an almost exponential behavior at high pj, and a flattening of the spectra 
at low pj. These features are compatible with the development of a strong collective flow with 
centrality, which dominates the spectral shapes up to relatively high pi in central collisions. The 
(j&r) parameter extracted from fits to the blast-wave parametrization indicates a radial flow about 
10% higher than at RHIC at a/^nn = 200 GeV in central collisions. The integrated abundances 
of particles are almost independent of centrality. They are compared with expectations from 
thermal models. While the K/7T ratio was found to agree with these expectations, the p/n 
is a factor 1.5 lower. The central collision data are successfully described by hydrodynamic 
models, but the low proton yield requires a refined description of the late fireball stages. These 
models ll47ll48Tl indicate a non-negligible baryon annihilation in the hadronic phase, which alters 
the thermal yields and leads to a lower proton yield. The origin of this effect, however, is not 
yet established and alternative explanations, such as non-equilibrium effects or flavor-dependent 
freeze-out, exist in the literature. In more peripheral collisions, purely hydrodynamic models 
give a poor description of the data, indicating the limit of applicability of hydrodynamics. 
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